Diproton correlation in a proton-rich Borromean nucleus 7 Ne 
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We carry out three-body-model calculations for a proton-rich Borromean nucleus 17 Ne by assum- 
ing a 15 + p + p structure. To this end, we use a density-dependent contact interaction between 
the valence protons, explicitly treating also the Coulomb interaction. We find that the two-particle 
density distribution for 17 Ne is similar to that for 16 C, which has two valence neutrons outside the 
N = 8 core. That is, the two protons take a spatially compact configuration, while the Coulomb 
repulsion plays a minor role. This indicates that there is a strong diproton correlation in the ground 
state of the 1 Ne nucleus. We also show that the Coulomb interaction reduces the expectation value 
of the proton-proton interaction by about 14%. 



PACS numbers: 21.10.Gv,21.60.-n, 21.45.+V, 27.20.+n 
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I. INTRODUCTION 

Physics of unstable nuclei have attracted much atten- 
tion both experimentally and theoretically [l|-y]- These 
nuclei are located far from the /3-stability line, having a 
large asymmetry in the proton-to-neutron ratio. In the 
neutron-rich side, a few exotic features have been ob- 
served. These include i) a large extension of density dis- 
tribution, referred to as halo or skin structure [4j, ii) a nar- 
row momentum distribution [5j, and iii) a large concentra- 
tion of the dipole strength distribution at low energies [g- 
8]. For n Li and 6 He nuclei, the so called Borromean 
structure has also been discussed extensively [9J, llOj . The 
Borromean is defined as a three-body bound system in 
which any two-body subsystem does not bound. The 
pairing interaction between the valence neutrons plays 
an essential role in stabilizing these nuclei [9J]. 

One of the topics of current interest in physics of unsta- 
ble nuclei is dineutron correlations in neutron-rich nuclei. 
Although two neutrons do not form a bound state in vac- 
uum, theoretical calculations have shown that they take a 
spatially compact configuration in finite nuclei [3|, iQl— ll 5jj . 
This feature is enhanced significantly when the binding 
is weak as in 11 Li and 6 He nuclei. The recent experimen- 
tal observation of the strong low-lying dipole strength 
distribution in the 11 Li nucleus [8J has provided an ex- 
perimental signature of the existence of strong dineutron 
correlation in this nucleus. 

In this paper, we discuss whether there exists a sim- 
ilar correlation in proton-rich nuclei. In contrast to the 
neutron-rich nuclei, proton-rich nuclei have been investi- 
gated less extensively, and the existence of diproton cor- 
relation has not yet been fully clarified. For this purpose, 
we choose 17 Ne nucleus. This nucleus is known to be a 
candidate of two-proton halo nucleus, as has been sug- 
gested by measured large interaction cross sections [16| as 
well as a recent measurement of charge radius [lTj . This 
nucleus is also known to have a Borromean nature, as its 
subsystem 16 F is proton unbound. The two-proton emis- 
sion decay has been observed from the first excited state 
of 17 Ne nucleus [18H20J . which suggests that the structure 



of 17 Ne can be understood as a three-body system of 
15 0+p + p. Based on this idea, several three-bodymodel 
calculations have been performed for 17 Ne [2lM26|. In 
this paper, we shall carry out similar three-body model 
calculations as in Refs. [11] . in which the dineutron cor- 
relation in n Li and 6 He has been investigated. This is 
based on a three-body model with a density-dependent 
pairing interaction between the valence nucleons [9|, [27|| . 
For the 17 Ne nucleus, we also include the Coulomb in- 
teraction explicitly and discuss the role of the Coulomb 
repulsion in diproton correlation. 

The paper is organized as follows. In Sec. II, we set up 
the three-body model with a density-dependent contact 
interaction. In Sec. Ill, we apply the model to the ground 
state of 17 Ne. We study the density distribution of the 
valence protons and discuss the diproton correlation. We 
also compare our results with those for the 16 C nucleus, 
whose neutron structure (N — 10) may have some simi- 
larity to the proton structure (Z = 10) in 17 Ne. We then 
summarize the paper in Sec. IV. 



II. THREE-BODY MODEL 

We describe the ground state of 17 Ne nucleus by as- 
suming an inert core nucleus 15 and two valence pro- 
tons. We consider the following three-body Hamiltonian 
for this system: 

tf(n,r 2 ) - / l W+/ l (2) + PilEi + ]/ pp ( ri)r2 ) ) (i) 

J\c r fTt 



^ = t+^), 



(2) 



where m and A c are the nucleon mass and the mass num- 
ber of the inert core nucleus, respectively, h is the single- 
particle (s.p.) Hamiltonian for a valence proton, in which 
V pc is the potential between the proton and the core nu- 
cleus. Vp P is the interaction between the valence protons. 
The diagonal component of the recoil kinetic energy of 
the core nucleus is included in the s.p. Hamiltonian h 
through the reduced mass \i = A c m/(A C + 1), whereas 
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FIG. 1: A low-lying level scheme of 16 F nucleus. Spin aver- 
aged levels for two sets of resonance states are also shown. 



the off-dia gona l part is taken into account in the last term 
inEq. g)|li|. 

For the single-particle potential Vp C , we use a Woods- 
Saxon potential together with the Coulomb interaction 
between a valence proton and the core nucleus. That is, 



V pc (r) 



r ar 



with 



m = 



exp [(r - -Score) /(We)] ' 



/(r) + V c (r), (3) 



(4) 
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where i? C ore = foA c is the radius of the core nucleus 
and Vq is the Coulomb potential for a uniform charge 
distribution. For the parameters, we use ro=1.22 fm, 
a corc =0.65 fm, V = -53.73 MeV, and V ls = 15.01 MeV. 
This parameter set yields an sjya and a d 5 / 2 resonances 
at 0.675 MeV and 1.129 MeV, respectively. These are 
spin-averaged energies of the observed _ resonance at 
0.535 MeV and 1~ resonance at 0.722 MeV, and 2" at 
0.951 MeV and 3" at 1.257 MeV in 16 F, respectively [l9| 
(see Fig. 1). 

We solve the three-body Hamiltonian, Eq. (JTJ, by ex- 
panding the ground state wave function with the single- 
particle wave functions, (i>rWjrn(r)[9.llll,l27j. 
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where we have assumed that the valence protons form 
a + configuration. The continuum s.p. states are dis- 
cretized in a large box with the size of -R D ox=30 fm. We 
include the s.p. angular momentum I up to 12. In the 
expansion, we explicitly exclude the lsi/2, l.P3/2j arL d 



l.Pi/2 states, which are occupied by the protons in the 
core nucleus. 

In this work, we employ a density-dependent contact 
interaction for the nuclear part of the proton-proton in- 
teraction [3, LkH [23, |30] , while the Coulomb potential is 
also explicitly included. That is, 



V PP (ri,r 2 ) = KW(n,r 2 ) + V^n,^), 



(7) 



with 



V p W(r 1 ,r 2 )=<5(r 1 -r 2 ) ,•„ 

and 

VgHr u r 2 ) = 
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As is well known, a contact interaction has to be han- 
dled in a truncated space defined with a cutoff energy 
E C \it- We take E cut = 60 MeV, so that the s.p. space is 
truncated as 12711 
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where e n ij is a s.p. energy. For the Coulomb interaction, 
we also use the same cutoff energy. We have confirmed 
that our results do not change significantly even if we use 
a larger value of E cut for the Coulomb interaction. For a 
given cutoff, the strength of the pairing interaction, Vo, 
is determined so as to reproduce an empirical scattering 
length a[27J], that is, 
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where k cut = ^/mE cut /h 2 . The empirical scattering 
length for a neutron- neutron scattering is a nn = —18.5 
fm [31( , while that for a proton-proton scattering is a pp — 
—7.8063 fm [321 ] . The difference between a nn and a pp is 
mainly due to the Coulomb repulsion in a two-proton sys- 
tem. Since we explicitly include the Coulomb interaction 
in our calculations, we use the neutron-neutron scatter- 
ing length a nn to determine the strength of the pairing 
interaction, Eq. (fTT|) . assuming charge independence of 
nuclear force. Once vq is determined in this way, the pa- 
rameters for the density-dependent term, v p ,R p , and a p 
are adjusted so as to reproduce the ground state energy 
of 17 Ne, -Eg.s. = —0.944 MeV, measured from the thresh- 
old of two-proton emission. See the first row in Table I 
for the values of the parameters. 



III. RESULTS 

We now numerically solve the three-body Hamiltonian 
and discuss the ground state properties of the 17 Ne nu- 
cleus. The results of our calculation are summarized in 



Table HIl As we mentioned in the previous section, these 
results were obtained with a = a nn = —18.5 fin. On 
the other hand, Refs. pH, [27J used a somewhat different 
value, a = —15 fm, in order to take into account finite 
momentum effects. We have checked that the results are 
almost the same even if we use a = — 15 fm, together 
with the parameter set for the pairing interaction shown 
in Table IT] 

In Table HH (Vp P ) and (Vp P ) are expectation val- 
ues of the nuclear and the Coulomb interactions, re- 
spectively. We find that the ratio of these quantities 
is about —0.14, that is, the Coulomb repulsion reduces 
the pairing energy by about 14%. This value is simi- 
lar to what has been found [28] using a non-empirical 
pairing energy density functional for proton pairing gaps 
29]. P([s 1/2 ] 2 ), P([d 5/2 ] 2 ), and P([d 3/2 } 2 ) are the prob- 
abilities of the Si/2,^5/2: and d 3 / 2 configurations in the 
ground state wave function, respectively. One sees that 
the d 5 / 2 configuration dominates in the ground state 



wave function. V^nn) = y / {ty g . s .\(r 1 



v^ 



r 2 ) 2 |^ g . s .) and 

_ 2N ) = v / (*g-s-l( r i+ r 2)74|*g.s.) are the root- 
mean-square (rms) distances between the valence protons 
and that between the center of mass of the protons and 
the core nucleus, respectively. With these rms distances, 
the increment of the radius of 17 Ne due to the valence 
protons is calculated as [23, [3(| , 



S(r 2 ) = (r 2 ) 



A c 



A -+ 2 A c + 2 
2A 



(r 2 ). 
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Once S (r 2 \ is evaluated in this way, the matter radius 
for 17 Ne can be estimated with Eq. (fl2|) to be 2.62 fm. 
To this end, we used an empirical rms matter radius of 
the core nucleus, r rms ( 15 0) =2.44 fm [38|. The s-wave 
probability, -P([si/2] 2 ) ; in our calculation is somewhat 
smaller than that in other calculations, e.g., Ref. [23J . 
and the rms radius for 17 Ne is slightly underestimated as 
compared to the experimental value, 2.75 fm 38]. 

(012 ) in Table II is the opening angle between the va- 
lence protons. The value of (#12) = 76.64 degrees is in a 



TABLE I: Parameter sets for the pairing interaction defined 
by Eq.(J5| used for the ground state of 17 Ne nucleus, a is the 
scattering length which determines the strength of the pairing 
interaction according to Eq. (1 1 If) . The cutoff energy is set to 
be -E cut =60 MeV. Each parameter set yields the two-proton 
separation energy of S2 P = 0.944 MeV. 





Vo (MeV) 


v p R p (fm) 


a p (fm) 


a= -18.5 fm 




-617.46 


-v 3.0912 


0.67 


a = -15.0 fm 




-608.37 


-Vo 3.0653 


0.67 


a = -7.81 fm, v£p } 


= 


-567.73 


-vo 3.1495 


0.67 



good agreement with that in Ref. [33[. In Ref. [33|, the 
opening angle was estimated using the total B(E1) value. 
We therefore compute also the sum rule value for the El 
transition (see Table I), 



B. 



SR (E1) = -Z 2 s e 2 <r 2 N _ 
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where 



Z c e = 



A c 



A c 



(14) 



(15) 



for a core+2p system, Z c being the proton number of 
the core nucleus (notice that Z c s = Z C /(A C + 2) for a 
corc+2n system). This factor is Z e ff = 0.41 for 17 Ne. 
We remind the readers that the actual value for the total 
B(E1) strength is somewhat smaller than the sum rule 
value due to the Pauli forbidden transitions [34| . 

It is instructive to compare our results with those ob- 
tained by mocking up the Coulomb interaction effectively 
in the density dependent contact interaction, i.e., by 
readjusting the parameters for the contact interaction 
without treating the Coulomb interaction explicitly. For 
this purpose, we use the proton-proton scattering length, 
a pp = —7.81 fm, to determine the strength of the pair- 
ing interaction vq. The other parameters are adjusted to 
reproduce the ground state energy. See the last row in 
Table [I] for the values of the parameters. The results are 
summarized in the second column in Table |TTJ One sees 
that the results are almost the same between the first 
and the second columns in Table [ill This indicates that 
the effect of Coulomb interaction on the pairing proper- 
ties can be well simulated by adjusting the parameters 
for the pairing interaction. 

We also compare our results for 17 Ne with those for 
16 C. Both nuclei have two valence nucleons outside the 
N or Z = 8 core, and one may expect some similari- 



TABLE II: Ground state properties for 17 Ne and 16 C ob- 
tained with a three-body model. See the text for the defini- 
tion for each quantity. 





17 Ne 


17 Ne (no Coul.) 


i6 C 


(V P < P ">) (MeV) 


-3.26 


-2.76 


-3.88 


(V^ } ) (MeV) 


0.448 








P([si/2] 2 ) (%) 


15.16 


15.91 


20.69 


P([d 5 /2\ 2 ) (%) 


75.19 


75.68 


64.97 


P([d 3 /2] 2 ) (%) 


3.83 


3.26 


6.63 


(r% N ) 1/2 (fm) 


4.688 


4.749 


4.579 


(r 2 c-2 N ) 1/2 (fm) 


3.037 


3.037 


3.099 


«S(r 2 ) 1/2 (fm) 


1.267 


1.273 


1.306 


(O12} (deg) 


76.64 


76.03 


74.35 


B(E1)sr (e 2 fm 2 ) 


1.494 


1.493 


1.290 


r( A Z) (fm) 


2.62 


2.62 


2.52 



ties between the two nuclei. The three-body model cal- 
culations for 16 C have been already performed in Refs. 

35J-437T] - We here repeat the same calculation as in Ref. 

35( , but using the density-dependent contact interaction. 
The results for 16 C are summarized in the last column 
in Table HI] In order to calculate the matter radius, we 
use r rms = 2.30 fm [3_8(| for the radius of the core nucleus 
14 C. Interestingly, the ground state properties are almost 
the same between 17 Ne and 16 C, despite that the ground 
state energy for 16 C (E„ s = —5.47 MeV) is significantly 
different from that for "Ne (E g . s , = -0.944 MeV). Both 
are d-wave dominant, and the opening angle for the va- 
lence nucleons is about 75 degrees. This similarity may 
be naturally understood if the the Coulomb interaction 
between the valence protons and the core nucleus in 17 Ne 
mainly perturbs the ground state energy without affect- 
ing significantly the ground state wave function. 

Let us next discuss the density distribution of the va- 
lence nucleons, 



p(ri,r 2 ,0 12 ) =\ * g . s .(ri,r 2 ,6»i 2 ) | 5 



(16) 



which arc normalized as 



/>00 />OC pTV 

/ 4nrldri / r\dr\ I 27rsin0i 2 d0i2 
Jo Jo Jo 

x P (r 1 ,r 2 ,8 12 ) = 1.(17) 

See Refs. [9|, [ll| for its explicit form. Figure 2 shows 
the density distribution for 17 Ne (the upper panel) and 
16 C (the lower panel). For a presentation purpose, we 
set the radius of the valence nucleons to be the same 
(i.e., r\ = r 2 = r) and multiply a weight factor of 
4-7rr 2 ■ 27rr 2 sin# 12 . One clearly sees that the density dis- 
tributions are similar to each other for both nuclei. That 
is, the density distribution is largely concentrated at a 
small opening angle dn ~ 20 (deg) (notice that the den- 
sity looks zero in the figure at Qyx — because of the 
weight factor). This implies that the Coulomb repulsion 
between the valence protons plays a minor role in the 
density distribution, causing a strong diproton correla- 
tion in the ground state of 17 Ne, as has been conjectured 
in Ref. [39| using a quasi-2D system. 

The three-peaked structure in the density distribution 
is due to the [ds/ 2 ] 2 configuration, that dominates in the 
ground state wave function. Figure 3 shows the density 
distribution for the pure [ds/ 2 ] 2 configuration for 17 Ne. 
This density distribution has symmetric three peaks, 
which is given by p(0) ex (5/4 • cos 4 9-1/2- cos 2 6 + 3/20) 
[35l | (the formula in Ref. [35| contains a small error, which 
is corrected here). With the pairing interaction, several 
single-particle levels with opposite parity are mixed up in 
the wave function, leading to the asymmetric distribution 
shown in Fig. 2. 

We demonstrate the strong diproton correlation in 
17 Ne also in another way. Figure 4 shows the density 
distribution of the second valence nucleon when the first 
nucleon is located along the z axis at the mean separa- 
tion distance between the core and the center of mass of 
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FIG. 2: (Color online) The density distribution for the ground 
state of Ne (the upper panel) and of C (the lower panel). 
It is plotted as a function of r = n = r2 and the opening 
angle between the valence nucleons, #12 • A weight factor of 
4nr 2 ■ 2nr 2 sin #12 has been multiplied. 
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FIG. 3: (Color online) The same as the upper panel of Fig. 
2, but for the pure [rfs/2] 2 configuration. 




the two valence nucleons, that is, r\ = \/( r c-2N)- ^ 
is plotted as a function of the coordinate of the second 
nucleon, z 2 = r 2 cosf?i 2 and x 2 = r 2 sin(9i 2 . As has been 
shown in Fig. 2, the second nucleon is well localized in 
the vicinity of the first nucleon, both for 17 Ne and 16 C 
nuclei. Evidently, the strong diproton correlation exists 
in the ground state of 17 Ne in spite of the Coulomb re- 
pulsion between the valence protons. 
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The upper and the lower panels are for 17 Ne and 16 C, respec- 
tively. 



this end, we used a density-dependent contact pairing in- 
teraction between the valence protons. In addition, we 
explicitly treated the Coulomb interaction as well. We 
found that the Coulomb repulsion leads to an about 14% 
reduction of the pairing energy. The density distribu- 
tion for the valence protons resembles the two-neutron 
density distribution in 16 C, showing a clear and strong 
diproton correlation in the ground state of 17 Ne. That is, 
the two valence protons take a spatially compact config- 
uration despite of the Coulomb repulsion. This suggests 
that the Coulomb interaction plays a minor role in the 
pairing correlation. In fact, we found that the effect of 
the Coulomb interaction can be well simulated by read- 
justing the parameters for the pairing interaction. 

The present study suggests that there is a strong dipro- 
ton correlation in light proton-rich nuclei. It would be 
an interesting future work to study whether this is the 
case also for medium- heavy and heavy proton-rich nuclei. 
The effect of the Coulomb interaction is expected to be 
stronger in the mean-field potential for these nuclei, but 
it is an open question whether the Coulomb effect on the 
diproton correlation is stronger or not. Another future 
work would be to study the role of diproton correlation 
in two-proton radioactivities. We will report on it in a 
separate paper. 
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